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Abstract 

In a two flavour world, the usual equal p, equal E, or equal v assump- 
tions used to derive the neutrino oscillation length can be replaced by a 
specific definition of the velocity of the oscillating state which, unlike those 
assumptions, is compatible with exact production kinematics. This defi- 
nition is further vindicated by the analysis of ir^ decay at rest. 
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1 Introduction 



Contradictory arguments are still exchanged in the discussion of neutrino oscil- 
lations in vaccuum not to mention more complicated cases. People keep arguing 
about the superiority of their favourite assumption (equal p's, E's, or u's ||) in 
deriving an oscillation length. Some of the champions of the last approximation 
even claim that the usual formula must be corrected by a factor of two Q while 
defenders of the standard result have spent some time showing that the equal 
v's is kinematically untenable . 

Although this last assumption is more at variance with exact production 
kinematics than the other two, treating it with care helps to uncover a weaker 
hypothesis which is both sufficient and compatible with kinematics, at least in 
a two flavour - two mass neutrino world. 

Using the naive plane-wave model in this simple case, we shall show that: 

1. The equal v assumption does yield the same oscillation length as the other 
two provided velocity is treated consistently, i.e. the same value is used 
in the momentum-energy and in the position-time relations. 

2. To derive this oscillation length, the only necessary assumption which is 
compatible with exact kinematics, and is, in some way, contained in the 
three usual approximations is a specific definition of the group velocity of 
the would-be wave packet. 

3. This value of the group velocity must be postulated, but is vindicated 
by the equality of the derived oscillation length with that found by more 
elaborate methods; moreover, it seems difficult to avoid in the case of 
7r — > fiu at rest. 

2 Standard and 'equal velocities' assumptions 

Being confined to a two mass world, we shall use simplified notations from the 
start and represent the neutrino born at space-time point (0, 0) in some charged 
current reaction involving charged lepton I = e or [i by: 

|0,0 >= \vi >= cob 6> 1 1 > + sin 6" 1 2 > (1) 

where \h > h = 1,2 are mass eigenstates with eigenvalues and definite 
energy and momentum Eh and ph- 

Assuming only the existence of a momentum operator commuting with the 
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Hamiltonian, (|T|) can be propagated by use of the space-time translation oper- 
ator U = e - l ( Ht - p ,F ) which yields: 

U\0,0 >= \x,t >= cosee' i( - Elt - pix) \l > +sme e - i ( E2t - p2x) \2 > (2) 

We have taken the x axis along the direction of motion. The values of Eh 
and ph depend on the production kinematics; in the case of a 7172 decay for 
example, they are fixed by the masses in the ir rest frame and from there in the 
lab, once the 7r's decay angle and velocity are given. Projection of \x,t > onto 
\vy >= — sin#|l > +cos#|2 > gives then 

Ai^ v {x,t) = sm9 cos 9(e-^ E2t - p2X ^ - e -^t~ Pl x)^ ^) 

for the amplitude to detect a neutrino of flavor V ^ I at point x and time t 
relative to the production point at (0,0), if that neutrino undergoes a C.C. 
interaction. The object of interest is the phase difference 5(f) = 5Et — 5px (with 
5E = E\ — E2 and 5p = pi — P2) which appears in the oscillation probability 
obtained by squaring (||): 

P^HM) = sin 2 (2#)sin 2 (^) (4) 

To evaluate 5(p, people postulate equal momenta, or equal energies, (see e.g. 
II) or even equal velocities (e.g. Q) of the two definite mass components. We 
first show that all three assumptions yield the same expression without further 
approximations : 



2.1 Equal energies 

It is assumed here that E\ = E2, hence 5<j) = —5px. However the relation: 

Sp 2 = 5E 2 - 5m 2 (5) 

yields in this case: 

5m 2 5m 2 5m 2 

op = = — and 0(f> = — —x 

pi + p 2 2p 2p 

with an obvious definition for p. 

Requiring then 5<j) = 2tt we find the well known result: 

Losc ~ \5ri r \ (6) 
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2.2 Equal momenta 



Here, 6(f) = SEt and (|5|) yields 6E = J^ E ■ If v is the velocity of the center of 
the would-be wave packet, one finds: 

dm 2 ft™ 2 /«N 

5$ = SEt = — — t = x (7 

V E x + E 2 2p v ; 

where the last equality holds exactly provided v = E ^ E 
With this definition for v, the same oscillation length as in (j2.l| ) is found without 
the traditional ultra-relativistic ingredients used to transform (Q) by setting 
t ~ x for the neutrino trajectory and E\ + E 2 ~ 2p . 



2.3 Equal velocities 



In this case, 5(f> does not reduce to a single term and it is very important not 
to "approximate" t by x, for in so doing one would arrive at: 

8(f) = (SE — 5p)x = 5me~ v x 

upon introducing r\ = tanh" 1 (t;) 



The oscillating pattern would be described by: 

I (^urii I 27TC^ 

sin 2 ( e~ v x) and the oscillation length : L' = — — - (8) 

2 \om\ 

However, since 

eV = Ex+pi = E 2 +p 2 _ 2{E + p) 
mi 77i2 mi + m2 

this yields finally 

_ 4(E+p)ir 

L o S c. - | 5m2 | iyj 

viz. twice the usual value in the relativistic regime where E ~ p; this is the 
origin of the claims of some authors. 

However, if the hypothesis of equal velocities has any meaning, then the cen- 
ter of the wave packet moves with that velocity, not with velocity 1. Therefore, 
defining its position by x = vt yields: 

SfTl 

6(f) = 5E(t — vx) = Srwyil/v — v)x = x (10) 

p| Now 

1 mi m-2 777-1 + "12 

777 pi p 2 2p 



1 7 = 1/vi — v 2 as usual 
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Hence 5<fi = ^r=-x and the formula found in results. 

Clearly, replacing 1/v — v by 1 — v (equivalent to t — ► x in |l^) cannot be 
harmless; in the situation of |2.2| , v is only an overall factor, and replacing it 
by 1 would induce a relative error on the phase difference which goes to with 
1 — v. Not so in the present case where the relative error is 1/(1 + u) - hence 
the factor 2 found above. Stated differently, the same value for v must be used 
in configuration and in momentum space. 

One sees that in || for example, the trap is escaped because of the further equal 
E's or equal p's assumption beyond the first t ~ x approximation. 
Note that even if the 'equal u's ' assumption is contradictory with production 
kinematics, as shown in [jlj] and 0, it must embody some approximate truth 
because, for oscillations to be observed, it is important that the two component 
waves do not separate too early, and therefore, that they do not have too 
different velocities. E\ sa E2 and, say, m\ <C m-2 as assumed in Q is possible 
only in ultra relativistic situations; but then very small velocity differences 
make huge differences in the 7 factor, so ^ <C 1 is possible, though exact 
equality is excluded; this last impossibility makes the loss of coherence and the 
disappearance of oscillations inescapable after a certain distance. 

3 A 'kinematically correct' assumption 
3.1 Requirements 

A connection between x and t is necessary to make predictions about observa- 
tions which will register a position w.r.t. the creation point, but never measure 
the elapsed time. This connection should arise from a full quantum treatment 
using wave packets or field theory; it would then only be necessary to propagate 
the system in time, with the space evolution being taken care of by the wave 
equation. 

In the simple plane wave approach, this relation must be imposed because plane 
waves are of infinite extent and with two different (mass) states involved one 
has apparently no mean of defining a 'center' either by averaging or by station- 
ary phase arguments. 

The three assumptions detailed above are just different ways of doing this. 
No further ingredient is needed in the equal v case, since there is an obvious 
definition for the group velocity. A very likely definition of this velocity has 
also been used in the equal p approximation, thereby avoiding the need of the 
'ultra relativistic' v ~ 1 and E\ + E2 ~ 2p . Finally, the equal E approximation 
circumvents the problem by yielding a stationary oscillation pattern from the 
beginning. ^ 

2 This makes it the preferred hypothesis of the authors of Q, [g| because it avoids ambigu- 
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It is clear now that if we endow the center of the would-be wave packet 
with the average momentum and energy of its components: p = P1 + P2 ) E = 
El + E2 and/or postulate that it has group velocity v g = we find back the 
same expression for the phase difference in all three cases; it appears that this 
postulate embodies all that is needed to yield the usual L osc : 

v = ^r^ 6<f> = 6p X -6Et = (6 P - SE^±^) X = ^ ~ 

Ei + E 2 pi + P2 Pi + P2 

—5m 2 x 



2p 



Indeed, this simple hypothesis is a/ implied by the ( |2.3|) assumption, b/ 
identical with the definition of the group velocity that we have used in the ( |2.2| ) 
case and c/ not required in the ( [2.1| ) derivation, but not incompatible either. 



3.2 The case of n — > [iv at rest 

We show here that our postulate is actually necessary in the kinematically most 
simple case of ir — ► /jlu at rest. 

Rrom Ej = - 2 ^ — and E„ + = m n one sees that: n 
mi < m 2 => Ei < E 2 =5- E^i > E^ 2 => P^i > P/j.,2 

and by momentum conservation pi > p 2 which, using again E\ < E 2 shows that 
v\ > v 2 ; as expected, all three assumptions are inconsistent with kinematics and 
the third is the worst since vi/v 2 = j^f^r 

Now, for oscillations to be observable, one must have \pi — p 2 \ < Ap where 
Ap is the quantal momentum uncertainty due to the localization of the source 
(the 7r in this case), so that one cannot tell which neutrino mass eigenstate 
has been produced. Consequently, the phase relation is preserved between the 
components of the muon state which can be described by a linear superposition 
of amplitudes. In the plane-wave approximation, this wave function must read 
something like: 

if, = cos6e- l{E ^ lt - pix) + sm6e~ l(E ^ 2t - p2x) 

where pi and p 2 are the same as for the two neutrino mass states recoiling 
against 11. 

ities in transforming theoretical results into observable spatial oscillation patterns. 
3 Which shows that 8E = contrary to the approximation of ^ 
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Therefore \ip\ 2 = 1 + sin 29 cos(5E fl t — 5px) and (assuming 9 < tt/2) maximizing 
the probability gives x = -j^-t = v^t for the center of the (very elementary) 
wave packet, with suggestive of a discrete version of the usual v g = 
However in the present case, = E Vl ^ p ^ - = 

But this makes the above postulate for the neutrino velocity more likely, because 
to find this last value for Vu, we have to endow the muon state with the mean 
momentum p and energy En of its two components, and likewise for the neutrino 
in order to enforce energy and momentum conservation in the decay; in turn, 
this implies our neutrino velocity. Of course, other ways of averaging could do 
the job as far as conservation laws are concerned; but it is hard to see how 
^-dependent averages of E and p could yield back the same (^-independent) 
expression for the velocity of the muon obtained by the requirement of maximum 
probability. ^ 

Certainly, assuming Ap > \pi — j?2 1 and using a superposition of two plane waves 
with p = pi and p = P2 for the muon is not perfectly consistent, but not worse 
either than what is done on the neutrino side. The reason for which we had to 
use the muon argument is evidently that the ill-defined oscillating state (||) has 
a constant modulus. Is is not possible to make any direct statement about its 
location in space-time without first defining its flavour content (as is done to 
calculate the oscillation length). 



4 Restrictions and conclusion 



The shortcomings of the plane wave description of neutrino oscillations were 



discussed long ago ; the necessity of using wave packets (see e.g.[|13|) or field 



theory ([14], [15]) and of incorporating the neutrino production and/or detection 



processes in the description of the phenomenon (|lC||, ||lq| ) has been the subject 
of many works and still feeds a continuing debate. However, the results of these 
more elaborate treatments always imply, up to damping factors, the oscillation 
pattern described by (||) in the relativistic limit, whenever oscillations are not 
washed out by resolution or loss of coherence []. 

The above analysis explains the robustness of this classical formula by ex- 
posing the only ingredient necessary to retrieve it; this turns out to be more 
general than what is used in standard derivations and does not contradict pro- 
duction kinematics. 

As a by-product, it also shows that the variable to be used in the oscillation 
length formula is the momentum rather than the energy, should the distinction 
apply. However this would happen in the non-relativistic regime where, due to 
phase-space limitations and admixture of 'wrong' helicity states, an oscillation 
probability disconnected from the rest of the process has a very limited meaning 

4 Which is but a simplified and discretized version of standard stationary phase arguments 
used to define the center of a wave packet 

5 Of course we discount incorrect treatments of the equal velocity assumption 
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0, see also @. 



Obviously, the restriction to two flavours is a weakness of the analysis pre- 
sented here; it must be remarked, however, that the oscillation length has a 
clear physical interpretation only in this case. The corresponding parameters 
become unavoidably entangled with mixing matrix elements when dealing with 
more than two flavours. Nevertheless, it is our hope that the simplicity of the 
arguments and results presented here may shed some light on the general case. 
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